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Abstract 

A numerical study was conducted to analyze the per- 
formance of different turbulence models when applied 
to the hypersonic NASA P8 inlet. Computational re- 
sults from the PARC2D code, which solves the full 
two-dimensional Reynolds-averaged Navier-Stokes equa- 
tions, have been compared with experimental data. The 
zero-equation models considered for the study were the 
Baldwin-Lomax model, the Thomas model, and a com- 
bination of the Baldwin-Lomax and Thomas models; the 
two-equation models considered were the Chien model, 
the Speziale model (both low Reynolds number), and 
the Launder and Spalding model (high Reynolds num- 
ber). The Thomas model performed best among the 
zero-equation models, and predicted good pressure distri- 
butions. The Chien and Speziale models compared very 
well with the experimental data, and performed better 
than the Thomas model near the walls. 

Nomenclature 

k turbulent kinetic energy 

M Mach number 

p static pressure 

P p pitot pressure 

Pt total pressure 

Re Reynolds number 

T t total temperature 

Ur friction velocity, \Jt w / p w 


’National Research Council-NASA Research Associate. Mem- 
ber AIAA 


X axial distance from the leading edge 

of the centerbody 

XREF inlet cowl height 

y + law-of-the-wall coordinate, u t Y/i> w 

Y vertical distance from the centerbody 

6 boundary-layer thickness 

e dissipation rate of turbulence energy 

v kinematic viscosity 

p density 

t w wall shear stress 


Subscripts: 

w evaluated at wall 

oo tunnel free stream condition 


Introduction 

Future hypersonic cruise vehicles will require a highly 
integrated and efficient propulsion system. Much of the 
success will depend on the progress made in computa- 
tional fluid dynamics for high-speed flows. Hypersonic 
flows are characterized by the presence of high Mach 
numbers in the inviscid region, high gradients of the flow 
variables in the vicinity of shocks, real gas effects, and 
thick heated boundary layers. To accurately model such 
complex flows, one needs to consider new turbulence and 
transition models, equilibrium and non-equilibrium real 
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gas effects, and the capability to resolve high flow gradi- 
ents. At present, a number of very powerful and efficient 
computational codes are available to solve the full Navier- 
Stokes equations. However, the success of these codes 
depends on, among other factors, the type of turbulence 
modeling being implemented in them. The objective of 
the present paper is to compare the performance of a 
variety of turbulence models for the prediction of hyper- 
sonic inlet flows. 

The NASA P8 inlet, 1 which represents cruise condi- 
tions of a typical hypersonic air-breathing vehicle, was se- 
lected as a test case for this study. Recently, an AGARD 
working group 2 also selected the P8 inlet for comparing 
numerical results. 

A relevant review of turbulence modeling is available 
in Ref. 3. In the past, Ng et al. 4 made an effort to 
compare the performance of different turbulence models 
for the calculation for a two-dimensional hypersonic in- 
let flow field and compared their results with experimen- 
tal data from the P8 inlet. They used the McDonald- 
Camarata model 5 and the Baldwin-Lomax 6 model (both 
zero-equation models) and the turbulence kinetic energy 
model 7 (one-equation model) in a parabolized Navier- 
Stokes (PNS) code, and obtained results for both fully 
turbulent flow throughout the inlet and while specify- 
ing the transition from laminar to turbulent flow. In 
their study, only the Baldwin-Lomax model was imple- 
mented in a full Navier-Stokes (NS) equation code. They 
concluded that the McDonald-Camarata model, with a 
user-specified transition of the boundary layer, gave the 
best solution among the three models. The Baldwin- 
Lomax model predicted separation for both the NS and 
PNS codes, where none existed in the experiments. How- 
ever, their study was not complete and conclusive since 
it was primarily limited to the PNS code, and only the 
Baldwin-Lomax turbulence model was implemented in a 
full Navier-Stokes equations code. Moreover, they did 
not include a two-equation turbulence model in their 
study. Therefore, there is a need to further investigate 
this problem in order to find the best existing turbulence 
model to predict the hypersonic inlet flow field. 

The PARC2D code, which solves the full two- 
dimensional Reynolds-averaged Navier-Stokes equations, 
was selected for this study. The PARC code has been 
used previously for a variety of supersonic and hypersonic 
configurations, 8,9 but only the Balwin-Lomax turbulence 
model was used. These investigations demonstrated the 
capability of the PARC code to simulate hypersonic inlet 
flows. 

The results are presented for six versions of zero- and 
two-equation turbulence models: the Baldwin-Lomax, 6 


Thomas, 10 and a combination of the Baldwin-Lomax and 
Thomas models, which are all zero-equation models; and 
the Chien, 11 Speziale et al. 12 (both low Reynolds num- 
ber), and Launder and Spalding 13 (high Reynolds num- 
ber) models, which are two-equation turbulence models. 

Experimental Background 

The experimental investigation was conducted at 
NASA Ames’ 3.5-Foot Hypersonic Wind Tunnel to de- 
termine the internal flow characteristics of a typical inlet 
on a hypersonic air-breathing vehicle operating at cruise 
conditions. The geometry of the inlet model tested is 
shown in Fig. 1. The inlet was a Mach 7.4 rectangular 
mixed-compression design with exiting supersonic flow. 
The model has an internal compression ratio of 8 and 
thus is referred to as the P8 inlet. Inlet cowl height 
(XREF) was 18.33 cm, forebody length was 82.28 cm, 
overall length was 136.2 cm, and width was 35.56 cm. A 
6.5° forebody wedge was designed to match an inlet en- 
trance Mach number of 6 at a test Mach number of 7.4, 
allowing for boundary layer displacement effects. The 
wedge was cooled to provide a uniform surface temper- 
ature of 0.375 T t00 , where the free stream toted temper- 
ature ( Tt oa ) was 8 11° K. The free stream unit Reynolds 
number (Re <») was 8.86 x 10 6 /m. The transition point 
on the centerbody was found experimentally to be at 
approximately 40 percent of the distance between the 
wedge leading edge and the inlet entrance. The transi- 
tion point in the cowl boundary layer was approximately 
halfway between the cowl leading edge and the throat 
station. The details of the experiments can be found in 
Ref. 1. The internal contours were designed to provide 
high total pressure recovery and uniform static pressure 
at the throat. To achieve this goal, it was necessary to 
cancel the cowl lip generated shock at the centerbody by 
appropriate design of the centerbody contour. The de- 
sign criteria of the shock cancellation at the centerbody 
was not satisfied experimentally because of inadequate 
treatment of the interaction of the cowl-generated shock 
and the turbulent centerbody boundary layer. The cen- 
terbody boundary layer was found to be of considerable 
size (approximately 15 percent of the inlet height at the 
point of the cowl shock and centerbody interaction). This 
interaction produced a downstream shock pattern that 
yielded significant nonuniformities at the inlet throat. 

This realistic inlet geometry, with strong viscous- 
inviscid interactions and the availability of extensive ex- 
perimental data, provided an excellent opportunity to 
verify the numerical algorithm and turbulence models 
for the application of predicting the hypersonic inlet flow 
field. 
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Turbulence Models 

Only a brief description of the turbulence models used 
for the present study is given here because of limited 
space. 

The Baldwin-Lomax model was chosen for the study 
because it is one of the most widely used turbulence mod- 
els. The model is an algebraic eddy viscosity turbulence 
model for separated flows, and its use avoids the neces- 
sity for finding the edge of the boundary layer. In this 
two-layer turbulence model, the vorticity generated near 
a solid wall is used to determine the length scales in the 
turbulent boundary layer. The details of the model are 
given in Ref. 6. 

Another turbulence model used in the present calcu- 
lations is also an algebraic eddy viscosity model and is 
loosely based on the Thomas formulation 10 of the Bald- 
win and Lomax model. 6 This model is referred to as the 
Thomas model in this paper. The details of the Thomas 
model can be found in Refs. 10 and 14. 

In the case of the combination of the Baldwin-Lomax 
and Thomas models, the code is designed to calculate 
the value of turbulent viscosity by both the Thomas and 
Baldwin-Lomax methods separately. Then it takes the 
larger of two values in the unbounded flow region and 
the smaller value near the no-slip boundary condition. 

In the present study, two low- Reynolds- number two- 
equation turbulence models were used. Since the wall 
shear stress is computed directly from the velocity gra- 
dient without having any preassumption of the velocity 
profile near the walls, the degree of empiricism is reduced 
considerably for these models. 

The Chien low-Reynolds-number model, 11 as modified 
by Nichols 15 for compressibility effects, was used in the 
present study. All the functions of Chien’s model are 
simple algebraic functions and do not involve any differ- 
ential operators. Therefore, this model is relatively easy 
to implement. 

The other low-Reynolds-number two-equation model 
used in this effort is by Speziale et al. 12 The numeri- 
cal solution for the k — e equations is that of Nichols. 15 
Modifications were made to improve the numerical sta- 
bility of the turbulence model and to extend the model 
to compressible flows. 

The high-Reynolds-number two-equation turbulence 
model used in the present study is by Launder and 
Spalding. 13 One of the assumptions made in deriving the 
wall functions for k and e is that the production and 
dissipation rates of the turbulence kinetic energy are in 


equilibrium. This assumption is strictly valid only in the 
logarithmic region of the boundary layer. Therefore, the 
first grid point must be placed in the logarithmic layer 
for correct use of the wall functions. The wall functions 
used in this model are essentially semi-empirical correla- 
tions connecting conditions at the first grid points and 
those at the wall. 

Method of Solution 

The PARC2D code, selected for this study, was de- 
veloped at the Arnolds Engineering Development Cen- 
ter. The PARC code, which is available in both two- 
and three-dimensional versions, solves the full Reynolds 
averaged Navier-Stokes equations in strong conservation 
form with the Beam and Warning approximate factor- 
ization. The code uses a central differencing scheme in 
a generalized curvilinear coordinate system with implicit 
and explicit second- and fourth-order artificial dissipa- 
tion. The details of the code can be found in Ref. 14. 

The computational grid used in this study was 221x91, 
nonuniform in the X direction. The grid was packed on 
both ends from the wedge leading egde to the cowl lead- 
ing edge, and was also geometrically stretched from the 
cowl leading edge to the outflow boundary. To resolve 
the viscous layer, the grid lines were packed close to the 
centerbody and cowl walls using hyperbolic tangent func- 
tions such that the first grid line was located, at a y + of 
approximately 1.0 away from the walls. As discussed 
earlier, the high-Reynolds-number turbulence model re- 
quires the first grid line to be in the logarithmic layer for 
proper use of wall functions. Therefore, a separate grid 
had to be made for the Launder and Spalding model. 
The grid was kept the same in the X direction, but a y + 
of approximately 30 was used away from the walls. The 
grid size was also kept the same (221 x 91) in this case. 

A nonreflective boundary condition, 9 using a simple 
Mach wave extrapolation, was applied on the upper 
boundary upstream of the cowl to let the shock wave 
from the wedge pass through the boundary. The flow 
field at the inflow boundary, which is ahead of the inlet 
wedge leading edge, was held fixed. The extrapolation 
boundary condition was applied to the outflow bound- 
ary. A no-slip isothermal boundary condition was used 
on all solid walls. 

It should be noted that all computations were per- 
formed with fully turbulent flow throughout the inlet. 
The transition points were not simulated in the present 
calculations because of code limitations. 

The computations were performed on the CRAY-YMP 
supercomputer at NASA Lewis Research Center. The 
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code typically required a total of 35 minutes of CPU time 
to achieve global convergence for the zero-equation mod- 
els. The computer time increased by approximately four 
times when the two-equation models were used. 

Results and Discussion 

The performance of a hypersonic inlet is significantly 
affected by the interaction of the shock waves and wall 
turbulent boundary layers. This interaction becomes 
more critical at higher Mach numbers, where boundary 
layer bleed is considered impractical because of the high 
stagnation temperatures. The interaction between the 
inlet cowl lip shock and centerbody turbulent boundary 
layer, particularly, requires careful analysis since the cen- 
terbody contours are often designed to cancel the cowl lip 
shock. The interned contours of the P8 inlet model were 
designed to provide cancellation of the cowl shock at the 
centerbody and an isentropic compression to the throat. 
Nevertheless, a reflected shock was found experimentally, 
which further interacted with the cowl and centerbody 
boundary layers upstream of the throat. This signif- 
icantly affected the flow structure at the inlet throat. 
The present computations were able to capture these flow 
characteristics. 

The computed density, pressure, and Mach number 
contours for the Thomas turbulence model are shown in 
Pig. 2. The contours for the other turbulence models 
were essentially similar to those of the Thomas model. It 
is clearly shown in the figure that the cowl shock, after in- 
teracting with the centerbody turbulent boundary layer, 
was reflected downstream and interacted with the cowl 
turbulent boundary layer. It was then reflected from the 
cowl and interacted again with the centerbody turbulent 
boundary layer before it left the inlet. 

The computational results are presented in the form 
of surface pressure distributions on the centerbody and 
cowl of the inlet model, and the pitot pressure and to- 
tal temperature distributions at many stations from the 
inlet entrance to the throat of the inlet. The results are 
compared with the corresponding experimented data of 
Ref. 1. 

The surface pressure distributions on the centerbody 
are shown in Fig. 3. The axial distances are non- 
dimensionalized with the inlet cowl height. As shown 
in the figure, the Thomas model compares very well 
with the experimental data, while the Baldwin-Lomax 
model and its combination with the Thomas model do 
not compare as accurately. The results obtained us- 
ing the Baldwin-Lomax model show the existence of a 
separation bubble on the centerbody in the immediate 


region of cowl-lip-shock/centerbody- boundary-layer in- 
teraction. No separation was reported in the experi- 
ment. Ng et al. 4 also reported the presence of a separa- 
tion bubble when they used the Baldwin-Lomax model. 
This indicates that the Baldwin-Lomax model fails in 
regions of shock- wave/boundary-layer interaction at hy- 
personic speeds. On the other hand, the Thomas, Chien, 
and Speziale models do not predict separation, but are 
able to successfully simulate the complex flow field re- 
sulting from the interaction of the cowl shock with the 
centerbody turbulent boundary layer. The Launder and 
Spalding model overpredicts the pressure rise at the inlet 
throat but predicts no separation on the centerbody. 

Figure 3 shows that all of the present computations 
disagree slighly with the experimented data near the lo- 
cation of the cowl shock wave and centerbody turbulent 
boundary layer interaction. This difference in the com- 
putational results with the experimental data may be at- 
tributed to the fact that the present computations were 
performed with fully turbulent flow throughout the in- 
let, without simulated transition points. As mentioned 
earlier, the transition point on the centerbody was found 
experimentally to be at approximately 40 percent of the 
distance between the wedge leading edge and the inlet en- 
trance. This may have resulted in thicker boundary lay- 
ers in the present computations as compared with those 
of experiments. The interaction of the cowl shock with a 
thicker turbulent boundary layer will produce a stronger 
upstream influence, which causes an upstream shift in 
pressure distributions. 

The pressure distributions on the cowl surface are 
presented in Fig. 4. The computed results show the 
flow compression and impingement of the reflected shock 
wave on the cowl. The expansion ahead of the re- 
flected shock impingement, which is a feature associated 
with the shock/boundary-layer interaction on the cen- 
terbody, is also visible in the numerical results. The 
Thomas model results are closest to the experimental 
data. The Baldwin-Lomax model and its combination 
with the Thomas model, and the Launder and Spald- 
ing model slightly overpredict cowl pressures at the inlet 
throat. However, the Chien and Speziale models slightly 
underpredict the cowl pressures at the throat. 

The computational and experimental pitot pressures 
at X/XREF = 5.67, a station upstream of the intersec- 
tion of the cowl shock with the centerbody, are shown 
in Fig. 5. The agreement of the numerical results with 
the experimental data is generally good. The steep rise 
in the pitot pressure at Y/XREF = 0.15 is due to the 
presence of the cowl shock. The computed results with 
the Thomas turbulence model seem to predict the best 
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solution among the zero-equation models. The Baldwin- 
Lomax model is not able to predict the inlet flow field 
accurately. All two-equation models are better than the 
zero-equation models near the centerbody. The Launder 
and Spalding model produces the sharpest cowl shock, 
and the calculated boundary layer thickness is much 
closer to that of experiment as compared with the other 
turbulence models. The results from the design analysis 
of Ref. 1 are also shown in the figure. The boundary layer 
thicknesses on the centerbody and cowl surfaces obtained 
from the experiments are also marked in the figure. 

All of the present computational results overpredicted 
the pitot pressures in the centred region of the inlet. In a 
recent AGARD report, 2 it was pointed out that no one 
has ever matched the experimental data using the stated 
tunnel conditions of Ref. 1. The many computational 
results tend to agree with one another, but do not match 
the experimented data. It was further suspected in Ref. 
2 that the conditions stated in Ref. 1 were different in 
some way from the conditions that were actually present 
in the tunnel. It should be noted that because of the 
high sensitivity of hypersonic flows, even small variations 
in the upstream flow field would lead to larger variations 
downstream. 

The total temperature profile at X/XREF = 5.67 is 
presented in Fig. 6. The experimental toted temperature 
ratio is less than 1.0 because it was not corrected for the 
uncertainties in other measurements 1 The comparison of 
the computed results with the experimental data is gen- 
erally good. The two-equation models compare very well 
with the experimental data, particularly near the center- 
body. 

The pitot pressure and toted temperature distributions 
at X/XREF = 6.09, the intersection point of cowl shock 
and centerbody, are shown in Figs. 7 and 8, respectively. 
The agreement between the computed results and exper- 
iment is good, except in the centred portion of the inlet, 
where the pitot pressures overpredict the experimented 
values. The performance of various turbulence models 
remains qualitatively the same as for X/XREF = 5.67. 

The pitot pressure and total temperature distributions 
at X/XREF = 6.37, a station just downstream of the re- 
flection of the cowl shock from the centerbody, are pre- 
sented in Figs. 9 and 10, respectively. As shown in Fig. 
9, the experimental centerbody boundary layer has been 
compressed by the reflecting shock wave and is compara- 
tively thinner than at the previous station. The reflected 
shock wave emerges from the boundary layer at Y/XREF 
= 0.025, as shown by the break in the curve. All of the 
present computations were able to pick up the emerging 
shock wave, but the magnitude and location of the emerg- 


ing shock varied for each turbulence model, as shown in 
Fig. 9. The design analysis of Ref. 1 is also presented 
in Fig. 9. The comparison of total temperatures with 
experimental data is good, as shown in Fig. 10. 

The pitot pressures and total temperatures at 
X/XREF = 6.65, a station upstream of the inlet throat, 
are shown in Figs. 11 and 12, respectively. The reflected 
shock from the centerbody interacts with the cowl bound- 
ary layer at X/XREF = 6.72, and gets reflected once 
again. Figure 11 shows that the centerbody boundary 
layer has been compressed thinner than at the previous 
station. As shown in Fig. 11, only the Thomas model was 
able to capture the presence of this reflected shock, while 
the other turbulence models including all two-equation 
models completely fail to do so. The comparison of to- 
tal temperatures with experimental data is qualitatively 
good, as shown in Fig. 12. 

The inlet performance was obtained in terms of pitot 
pressure, total temperature, and Mach number distri- 
butions at the throat. The pitot pressure distribution, 
presented in Fig. 13, shows a large variation across the 
throat height due to the presence of the reflected shock 
wave. The results from the design analysis of Ref. 1 are 
also shown in the figure. The total temperature distribu- 
tion at the throat is shown in Fig. 14. The Mach number 
profile at the throat is presented in Fig. 15. The com- 
puted results are in fair agreement with the experimented 
data. The results from the design analysis of Ref. 1 are 
also presented in the figure. 

Conclusions 

A computational study has been carried out to eval- 
uate the performance of various turbulence models in 
predicting hypersonic inlet flow fields. The PARC2D 
code, which solves the full two-dimensional Reynolds- 
averaged Navier-Stokes equations, has been selected to 
compare the computational results with available ex- 
perimental data. The results have been presented for 
the Baldwin-Lomax, Thomas, and a combination of the 
Baldwin-Lomax and Thomas models, which are all zero- 
equation models; and the Chien, Speziale (both low 
Reynolds number), and Launder and Spalding (high 
Reynolds number) models, which are two-equation tur- 
bulence models. 

These models have been evaluated on the basis of com- 
parisons of pitot pressure and total temperature profiles 
at various axial locations and surface pressure distribu- 
tions on the centerbody and cowl of the inlet. Based 
on the present results, the Thomas model compares very 
well with the experimental data, and it performs best 
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among the zero-equation models. The Baldwin-Lomax 
model and its combination with the Thomas model are 
not able to resolve the problem of shock wave and bound- 
ary layer interaction accurately. The Baldwin-Lomax 
model predicts separation near the interaction of the cowl 
shock with the centerbody boundary layer, where none is 
known to exist in experiments. The Chien and Speziale 
models compare very well with the experimental data, 
and perform better than the Thomas model, particu- 
larly near the walls. However, the Launder and Spalding 
model does not perform well as the Chien and Speziale 
models. Considering the fact that the CPU time required 
for the Thomas model is far less than the two-equation 
models, it is concluded that the Thomas model is best 
suited for the prediction of pressure distributions, and 
the Chien and Speziale (both-low-Reynolds number two- 
equation) models are recommended for the prediction of 
flow quantities near the walls and may be used for the 
calculation of skin friction and heat transfer coefficients. 
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Figure 1 . — Schematic diagram of P8 inlet model. 
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Figure 2. — Typical density, pressure, and Mach number contours for P8 inlet. 
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Figure 4. — Surface pressure distribution on the cowl. 
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Figure 9. — Pitot pressure distribution at X/XREF = 6.37. 
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Figure 10. — Total temperature distribution at X/XREF = 6.37. 
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Figure 11. — Pitot pressure distribution at X/XREF = 6.65. 
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Figure 12. — Total temperature distribution at X/XREF = 6.65. 
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Figure 1 3. — Pitot pressure distribution at the throat of 
P8 iniet. 
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Figure 14. — Total pressure distribution at the throat of 
P3 inlet. 
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Figure 1 5. — Mach number distribution at the throat of 
P8 iniet. 
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